Mechanism of DNA origami folding elucidated by mesoscopic simulations.

Many experimental and computational efforts have sought to understand DNA origami folding, but the time and length scales of this process pose significant challenges. Here, we present a mesoscopic model that uses a switchable force field to capture the behavior of single- and double-stranded DNA motifs and transitions between them, allowing us to simulate the folding of DNA origami up to several kilobases in size. Brownian dynamics simulations of small structures reveal a hierarchical folding process involving zipping into a partially folded precursor followed by crystallization into the final structure. We elucidate the effects of various design choices on folding order and kinetics. Larger structures are found to exhibit heterogeneous staple incorporation kinetics and frequent trapping in metastable states, as opposed to more accessible structures which exhibit first-order kinetics and virtually defect-free folding. This model opens an avenue to better understand and design DNA nanostructures for improved yield and folding performance.

Salt-Dependent Self-Association of Trinucleotide Repeat RNA Sequences.

Repeat RNA sequences self-associate to form condensates. Simulations of a coarse-grained single-interaction site model for (CAG)n (n = 30 and 31) show that the salt-dependent free energy gap, ΔGS, between the ground (perfect hairpin) and the excited state (slipped hairpin (SH) with one CAG overhang) of the monomer for (n even) is the primary factor that determines the rates and yield of self-assembly. For odd n, the free energy (GS) of the ground state, which is an SH, is used to predict the self-association kinetics. As the monovalent salt concentration, CS, increases, ΔGS and GS increase, which decreases the rates of dimer formation. In contrast, ΔGS for shuffled sequences, with the same length and sequence composition as (CAG)31, is larger, which suppresses their propensities to aggregate. Although demonstrated explicitly for (CAG) polymers, the finding of inverse correlation between the free energy gap and RNA aggregation is general.

Torsion affects the calculation of DNA twisting number.

The topological properties of DNA have long been a focal point in biophysics. In the 1970s, White proposed that the topology of closed DNA double helix follows White's formula: Lk=Wr+Tw. However, there has been controversy in the calculation of DNA twisting number, partly due to discrepancies in the definition of torsion in differential geometry. In this paper, we delved into a detailed study of torsion, revealing that the calculation of DNA twisting number should use the curve's geodesic torsion. Furthermore, we found that the discrepancy in DNA twisting numbers calculated using different torsion is N. This study elucidated the impact of torsion on the calculation of DNA twisting numbers, aiming to resolve controversies in the calculation of DNA topology and provided accurate computational methods and theoretical foundations for related research.

Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences.

Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

RNA: The Unsuspected Conductor in the Orchestra of Macromolecular Crowding.

This comprehensive Review delves into the chemical principles governing RNA-mediated crowding events, commonly referred to as granules or biological condensates. We explore the pivotal role played by RNA sequence, structure, and chemical modifications in these processes, uncovering their correlation with crowding phenomena under physiological conditions. Additionally, we investigate instances where crowding deviates from its intended function, leading to pathological consequences. By deepening our understanding of the delicate balance that governs molecular crowding driven by RNA and its implications for cellular homeostasis, we aim to shed light on this intriguing area of research. Our exploration extends to the methodologies employed to decipher the composition and structural intricacies of RNA granules, offering a comprehensive overview of the techniques used to characterize them, including relevant computational approaches. Through two detailed examples highlighting the significance of noncoding RNAs, NEAT1 and XIST, in the formation of phase-separated assemblies and their influence on the cellular landscape, we emphasize their crucial role in cellular organization and function. By elucidating the chemical underpinnings of RNA-mediated molecular crowding, investigating the role of modifications, structures, and composition of RNA granules, and exploring both physiological and aberrant phase separation phenomena, this Review provides a multifaceted understanding of the intriguing world of RNA-mediated biological condensates.

Application of Integrated Computational Approaches in Prediction of Plant Virus Encoded miRNAs and Their Targeted Plant Genes.

This chapter presents a comprehensive approach to predict novel miRNAs encoded by plant viruses and identify their target plant genes, through integration of various ab initio computational approaches. The predictive process begins with the analysis of plant viral sequences using the VMir Analyzer software. VMir Viewer software is then used to extract primary hairpins from these sequences. To distinguish real miRNA precursors from pseudo miRNA precursors, MiPred web-based software is employed. Verified real pre-miRNA sequences with a minimum free energy of < -20 Kcal/mol, are further analyzed using the RNAshapes software. Validation of predictions involves comparing them with available Expressed Sequence Tags (ESTs) from the relevant plant using BlastN. Short sequences with lengths ranging from 19 to 25 nucleotides and exhibiting <5 mismatches are prioritized for miRNA prediction. The precise locations of these short sequences within pre-miRNA structures generated using RNAshapes are meticulously identified, with a focus on those situated on the 5' and 3' arms of the structures, indicating potential miRNAs. Sequences within the arms of pre-miRNA structures are used to predict target sites within the ESTs of the specific plant, facilitated by psRNA Target software, revealing genes with potential regulatory roles in the plant. To confirm the outcome of target prediction, results are individually submitted to the RNAhybrid web-based software. For practical demonstration, this approach is applied to analyze African cassava mosaic virus (ACMV) and East African cassava mosaic virus-Uganda (EACMV-UG) viruses, as well as the ESTs of Jatropha and cassava.

Reversible strain-promoted DNA polymerization.

Molecular strain can be introduced to influence the outcome of chemical reactions. Once a thermodynamic product is formed, however, reversing the course of a strain-promoted reaction is challenging. Here, a reversible, strain-promoted polymerization in cyclic DNA is reported. The use of nonhybridizing, single-stranded spacers as short as a single nucleotide in length can promote DNA cyclization. Molecular strain is generated by duplexing the spacers, leading to ring opening and subsequent polymerization. Then, removal of the strain-generating duplexers triggers depolymerization and cyclic dimer recovery via enthalpy-driven cyclization and entropy-mediated ring contraction. This reversibility is retained even when a protein is conjugated to the DNA strands, and the architecture of the protein assemblies can be modulated between bivalent and polyvalent states. This work underscores the utility of using DNA not only as a programmable ligand for assembly but also as a route to access restorable bonds, thus providing a molecular basis for DNA-based materials with shape-memory, self-healing, and stimuli-responsive properties.

An RNA origami robot that traps and releases a fluorescent aptamer.

RNA nanotechnology aims to use RNA as a programmable material to create self-assembling nanodevices for application in medicine and synthetic biology. The main challenge is to develop advanced RNA robotic devices that both sense, compute, and actuate to obtain enhanced control over molecular processes. Here, we use the RNA origami method to prototype an RNA robotic device, named the "Traptamer," that mechanically traps the fluorescent aptamer, iSpinach. The Traptamer is shown to sense two RNA key strands, acts as a Boolean AND gate, and reversibly controls the fluorescence of the iSpinach aptamer. Cryo-electron microscopy of the closed Traptamer structure at 5.45-angstrom resolution reveals the mechanical mode of distortion of the iSpinach motif. Our study suggests a general approach to distorting RNA motifs and a path forward to build sophisticated RNA machines that through sensing, computing, and actuation modules can be used to precisely control RNA functionalities in cellular systems.

Coarse-grained modeling of DNA-RNA hybrids.

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

Pulley Effect in the Capture of DNA Translocation through Solid-State Nanopores.

Nanopores are powerful single-molecule sensors for analyzing biomolecules such as DNA and proteins. Understanding the dynamics of DNA capture and translocation through nanopores is essential for optimizing their performance. In this study, we examine the effects of applied voltage and pore diameter on current blockage, translocation time, collision, and capture location by translocating λ-DNA through 5.7 and 16 nm solid-state nanopores. Ionic current changes are used to infer DNA conformations during translocation. We find that translocation time increases with pore diameter, which can be attributed to the decrease of the stall force. Linear and exponential decreases of collision frequency with voltage are observed in the 16 and 5.7 nm pores, respectively, indicating a free energy barrier in the small pore. Moreover, the results reveal a voltage-dependent bias in the capture location toward the DNA ends, which is explained by a "pulley effect" deforming the DNA as it approaches the pore. This study provides insights into the physics governing DNA capture and translocation, which can be useful for promoting single-file translocation to enhance nanopore sensing.

Zα Domain of ADAR1 Binds to an A-Form-like Nucleic Acid Duplex with Low Micromolar Affinity.

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of viral and innate immune response proteins. While Z-form adoption is preferred by certain sequences, such as the commonly studied (CpG)n repeats, Zα has been reported to bind to a wide range of sequence contexts. Studying how Zα interacts with B-/A-form helices prior to their conversion to the Z-conformation is challenging as binding coincides with Z-form adoption. Here, we studied the binding of Zα fromHomo sapiens ADAR1 to a locked "A-type" version of the (CpG)3 construct (LNA (CpG)3) where the sugar pucker is locked into the C3'-endo/C2'-exo conformation, which prevents the duplex from adopting the alternating C2'/C3'-endo sugar puckers found in the Z-conformation. Using NMR and other biophysical techniques, we find that ZαADAR1 binds to the LNA (CpG)3 using a similar interface as for Z-form binding, with a dissociation constant (KD) of ∼4 µM. In contrast to Z-DNA/Z-RNA, where two ZαADAR1 bind to every 6 bp stretch, our data suggests that ZαADAR1 binds to multiple LNA molecules, indicating a completely different binding mode. Because ZαADAR1 binds relatively tightly to a non-Z-form model, its binding to B/A-form helices may need to be considered when experiments are carried out which attempt to identify the Z-form targets of Zα domains. The use of LNA constructs may be beneficial in experiments where negative controls for Z-form adoption are needed.

Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences.

Interplay between divalent cations (Mg2+ and Ca2+) and single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA), as well as stacking interactions, is important in nucleosome stability and phase separation in nucleic acids. Quantitative techniques accounting for ion-DNA interactions are needed to obtain insights into these and related problems. Toward this end, we created a sequence-dependent computational TIS-ION model that explicitly accounts for monovalent and divalent ions. Simulations of the rigid 24 base-pair (bp) dsDNA and flexible ssDNA sequences, dT30 and dA30, with varying amounts of the divalent cations show that the calculated excess number of ions around the dsDNA and ssDNA agree quantitatively with ion-counting experiments. Using an ensemble of all-atom structures generated from coarse-grained simulations, we calculated the small-angle X-ray scattering profiles, which are in excellent agreement with experiments. Although ion-counting experiments mask the differences between Mg2+ and Ca2+, we find that Mg2+ binds to the minor grooves and phosphate groups, whereas Ca2+ binds specifically to the minor groove. Both Mg2+ and Ca2+ exhibit a tendency to bind to the minor groove of DNA as opposed to the major groove. The dA30 conformations are dominated by stacking interactions, resulting in structures with considerable helical order. The near cancellation of the favorable stacking and unfavorable electrostatic interactions leads to dT30 populating an ensemble of heterogeneous conformations. The successful applications of the TIS-ION model are poised to confront many problems in DNA biophysics.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop.

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

High-Fidelity RNA Copying via 2',3'-Cyclic Phosphate Ligation.

Templated ligation offers an efficient approach to replicate long strands in an RNA world. The 2',3'-cyclic phosphate (>P) is a prebiotically available activation that also forms during RNA hydrolysis. Using gel electrophoresis and high-performance liquid chromatography, we found that the templated ligation of RNA with >P proceeds in simple low-salt aqueous solutions with 1 mM MgCl2 under alkaline pH ranging from 9 to 11 and temperatures from -20 to 25 °C. No additional catalysts were required. In contrast to previous reports, we found an increase in the number of canonical linkages to 50%. The reaction proceeds in a sequence-specific manner, with an experimentally determined ligation fidelity of 82% at the 3' end and 91% at the 5' end of the ligation site. With splinted oligomers, five ligations created a 96-mer strand, demonstrating a pathway for the ribozyme assembly. Due to the low salt requirements, the ligation conditions will be compatible with strand separation. Templated ligation mediated by 2',3'-cyclic phosphate in alkaline conditions therefore offers a performant replication and elongation reaction for RNA on early Earth.

Reverse transcription of plasma-derived HIV-1 RNA generates multiple artifacts through tRNA(Lys-3)-priming.

In vitro reverse transcription of full-length HIV-1 RNA extracted from the blood plasma of people living with HIV-1 remains challenging. Here, we describe the initiation of reverse transcription of plasma-derived viral RNA in the absence of an exogenous primer. Real-time PCR and Sanger sequencing were applied to identify the source and to monitor the outcome of this reaction. Results demonstrated that during purification of viral RNA from plasma, tRNA(Lys-3) is co-extracted in a complex with the viral RNA. In the presence of a reverse transcription enzyme, this tRNA(Lys-3) can induce reverse transcription, a reaction that is not confined to transcription of the 5' end of the viral RNA. A range of cDNA products is generated, most of them indicative for the occurrence of in vitro strand transfer events that involve translocation of cDNA from the 5' end to random positions on the viral RNA. This process results in the formation of cDNAs with large internal deletions. However, near full-length cDNA and cDNA with sequence patterns resembling multiple spliced HIV-1 RNA were also detected. Despite its potential to introduce significant bias in the interpretation of results across various applications, tRNA(Lys-3)-driven reverse transcription has been overlooked thus far. A more in-depth study of this tRNA-driven in vitro reaction may provide new insight into the complex process of in vivo HIV-1 replication.IMPORTANCEThe use of silica-based extraction methods for purifying HIV-1 RNA from viral particles is a common practice, but it involves co-extraction of human tRNA(Lys-3) due to the strong interactions between these molecules. This co-extraction becomes particularly significant when the extracted RNA is used in reverse transcription reactions, as the tRNA(Lys-3) then serves as a primer. Reverse transcription from tRNA(Lys-3) is not confined to cDNA synthesis of the 5' end of the viral RNA but extends across various regions of the viral genome through in vitro strand transfer events. Co-extraction of tRNA(Lys-3) has been overlooked thus far, despite its potential to introduce bias in downstream, reverse transcription-related applications. The observed events in the tRNA(Lys-3)-induced in vitro reverse transcription resemble in vivo replication processes. Therefore, these reactions may offer a unique model to better understand the replication dynamics of HIV-1.

A beginner's handbook to identify and characterize i-motif DNA.

Genomic DNA exhibits an innate ability to manifest diverse sequence-dependent secondary structures, serving crucial functions in gene regulation and cellular equilibrium. While extensive research has confirmed the formation of G-quadruplex structures by guanine-rich sequences in vitro and in cells, recent investigations have turned the quadruplex community's attention to the cytosine (C)-rich complementary strands that can adopt unique tetra-stranded conformation, termed as intercalated motif or i-motif. I-motifs are stabilized by hemi-protonated C:CH+ base pairs under acidic conditions. Initially, the in vivo occurrence of i-motifs was underestimated because their formation is favored at non-physiological pH. However, groundbreaking research utilizing the structure-specific iMab antibody and high-throughput sequencing have recently detected their conserved dispersion throughout the genome, challenging previous assumptions. Given the evolving nature of this research field, it becomes imperative to conduct independent in vitro experiments aimed at identifying potential i-motif formation in C-rich sequences and consolidating the findings to address the properties of i-motifs. This chapter serves as an introductory guide for the swift identification of novel i-motifs, where we present an experimental framework for investigating and characterizing i-motif sequences in vitro. In this chapter, we selected a synthetic oligonucleotide (C7T3) sequence and outlined appropriate methodologies for annealing the i-motif structure into suitable buffers. Then, we validated its formation by CD (Circular Dichroism) and NMR (Nuclear Magnetic Resonance) spectroscopy. Finally, we provided a thorough account of the step-by-step procedures to investigate the effect of i-motif formation on the stalling or retardation of DNA replication using high resolution primer extension assays.

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields.

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.

Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism.

G-quadruplexes (G4) are helical structures found in guanine-rich DNA or RNA sequences. Generally, their formalism is based on a few dozen structures, which can produce some inconsistencies or incompleteness. Using the website ASC-G4, we analyzed the structures of 333 intramolecular G4s, of all types, which allowed us to clarify some key concepts and present new information. To each of the eight distinguishable topologies corresponds a groove-width signature and a predominant glycosidic configuration (gc) pattern governed by the directions of the strands. The relative orientations of the stacking guanines within the strands, which we quantified and related to their vertical gc successions, determine the twist and tilt of the helices. The latter impact the minimum groove widths, which represent the space available for lateral ligand binding. The G4 four helices have similar twists, even when these twists are irregular, meaning that they have various angles along the strands. Despite its importance, the vertical gc succession has no strict one-to-one relationship with the topology, which explains the discrepancy between some topologies and their corresponding circular dichroism spectra. This study allowed us to introduce the new concept of platypus G4s, which are structures with properties corresponding to several topologies.